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Abstract 

We consider a stochastic perturbation of a FitzHugh-Nagumo system. 
We show that it is possible to generate oscillations for values of pa- 
rameters which do not allow oscillations for the deterministic system. 
We also study the appearance of a new equilibrium point and new 
bifurcation parameters due to the noisy component. 
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1 Introduction. 



Let us consider the following family of deterministic systems indexed by the 
parameters a G R and 6 > 0. 

6xt = -yt + f{xt), Xo = x (1.1) 
i/t = xt-a, Yo = y (1.2) 

and their stochastic perturbation by a one dimensional Wiener process {Wt) 
as follows 

6dXt = i-Yt + f{Xt))dt + y/~edWt, Xo = x (1.3) 
dYt = {Xt-a)dt, Yo = y (1.4) 

The function / is a cubic polynomial: f{x) = —x{x — a){x — (3) with 
a < < /3. The parameter 5 is small. The deterministic system f ll.ip -(1.2) 
is an example of a slow-fast system: the two variables x, y have different 
time scales, Xt evolves rapidly while yt evolves slowly. This system is one 
version of the so called FitzHugh-Nagumo system and plays an important 
role in neuronal modelling. In this context Xt denotes the voltage or ac- 
tion potential of the membrane of a single neuron. It was first proposed by 
FitzHugh and Nagumo (cf. [3], [I3]). One interest of this model is that it 
reproduces periodic oscillations observed experimentally. Indeed FitzHugh- 
Nagumo system finds its origin in the nonlinear oscillator model proposed by 
van der Pol. It is also a simplification of the Ho dgkin- Huxley model which 
describes the coupled evolution of the membrane potential and the different 
ionic currents: existence of different time scales enable to pass from a four 
dimensional model to a two dimensional one. Oscillations can take place be- 
cause the deterministic system fll.ip -(1.2) exhibits bifurcations; more details 
will be given in section 3. Let us mention that oscillations in this system 
f ll.ll) -(1.2) can only occur when a G]ao,ai[ where oq < ai are two particular 
values of parameter a namely the bifurcation parameters. 
Our main interest in the present paper is to generate oscillations even for 
a < ao (symmetrically a > ai) )by adding a stochastic perturbation to the 
deterministic system. What may be interpreted as some resonance type 
effect (cf. in], [H])- We will therefore investigate possible oscillations for 
system fll.3l) -(1.4). The presence of parameter e introduces a third scale in 
the system and the relative strength of 5 and e measured by the ratio ^i^-^i 
will determine its evolution. Our study was inspired by reference [6j where 
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M. Freidlin considers a random perturbation of the second order equation 
^li^ ~ ^^'^ performs the study of its solution using the theory of 

large deviations (cf. [5]). See also [7] for the study of a more general situa- 
tion. In our case g{y,y) = y — f{y + a). Although our argument is close to 
M. Freidlin's, the presence of parameter a leads to a richer behaviour. 
We prove the existence of equilibrium point and limit cycles different from the 
deterministic ones; as in [B] , as well as a new bifurcation point which did not 
exist for the deterministic system. Our study relies on transitions between 
basins of attraction of stable equilibrium points due to noise. Relying on 
some estimation of a family of exit times (propositions 3.3 and 3.4) we study 
conditions on the parameters under which a convenient stochastic dynamic 
system approach its main state (proposition 3.5 ), (which corresponds to the 
equilibrium point exhibited in main theorem 2.2), or approach a metastable 
state (which corresponds to the limit cycle and to the new bifurcation pa- 
rameters exhibited in main theorem 2.1). 

A general study of slow-fast systems perturbed by noise can be found in 
pp. Bursting oscillations in which a system alternates periodically between 
phases of quiescence and phases of repetitive spiking has been studied for 
stochastically perturbed systems in [10] and may be studied later in our 
stochastic setting. We recall that in the deterministic one a bursting-type 
behaviour has been generated in [2]. 

The paper is organized as follows. In section 2 we recall basic facts about 
f ll.ip -(1.2) and we state the two main theorems (2.1) and (2.2). Section 3 is 
devoted to the application of large deviation theory to fll.3p -(1.4) and section 
4 to the proof of the main theorems. 

2 Some Basic Results. 

2.1 Deterministic FitzHugh-Nagumo System. 

In f ll.ip -(1.2) let us consider a < < P and /(x) = —x{x — a){x — 13). In 
order to investigate the asymptotic behaviour of (xt, yt) one first looks for 
equilibrium points of the system and their stability. The equilibrium points 
are defined as the points (x, y) where the right-hand sides of both equations of 
the system vanish. For any value of a there is therefore a unique equilibrium 
point for fll.ip -(1.2) which is {a,f{a)). Moreover let ao,ai,ao < ai be the 
two points where /' vanishes. The stability of the equilibrium point (a, / (a)) 
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changes when a passes through the value Oq (resp. Oi); oq and ai are called 
the bifurcation parameters of the system. Let us focus on oq; an analogous 
argument holds for Oi. By linearizing system (II. ip at (ao + r], f{ao + r])) for 
1] small, we obtain system Z = AZ with 

/ /'(«o+r?) _1 \ 



A admits the two eigenvalues X± = ^(/'(cto + v) ^ — /'^(oq + rj)). 
The sign of /'(oq + 1]) is the same as that of rj since /' is increasing in the 
neighbourhood of Oq. The point (ag + 1], /(ag + is therefore an attracting 
(resp. repulsive) focus when rj < (resp. t] > 0). In particular, (a, /(a)) 
is stable when a < ao, unstable when a > gq. For a E]ao,ai[ the system 
admits a limit cycle. The bifurcation is of Hopf type [1] It can be verified 
numerically that if 6 < 0.01 the limit cycle is very close to the loop made 
up with the two attracting branches of the curve y = f{x) where x t— > f{x) 
is decreasing and y G [/(oq), /(ai)], and the portions of the two horizontal 
segments y = /(ao); V = /(^i) connecting them. When 5 — >• the period of 
this limit cycle is 0(1); by example for f{x) = a;(4 — x'^) it is equal to 2 (cf. 

)• 



2.2 Main Theorems 

Consider {Xf, Yf) the solution of (1.3)-(1.4) and S* > given in definition 3.3. 
Let us assume that e > and 5 > go to zero in such a way that for some 
constant c > 0. 

(2.1) 



Let us denote by lim * any limit on e and 6 going to under condition (12.1 p . 
Theorem 2.1 Let c g]0, S[; then: 

Ufa G]a:_(c), x+(c)[ where x_(c) and x^{c) are introduced in definition 3.4; 
then there exist two periodic functions and given in definition 4-2, s.t 
for all A, h>0, y e]/(ao), /(ai)[, 

lim* P(,,,) (^^|Xt-$«(t)|2dt>/i) = (2.2) 
lim* P(,,,) (sup \Yt - ^:(t)| >h) = (2.3) 

[0,A] 
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2. If a < a;_(c) or a > a;+(c) then for all y G]/(ao), for all h > 
there exists t{y,h) such that for all A > t{y,h), 

lim*P(^,^)( sup \Xt -a\ + \Yt - f{a)\ > h) ^ (2.4) 




Figure 1: Limit cycle when f (x) = x{4: — x'^) and — > c 
Remcirk 

In the first case the solution stabilizes when 5^0 and ^ii^^ c around 
a limit cycle defined by c and different from the one obtained in the deter- 
ministic case when 6^0 and e = (see figure 1). Moreover ,t_(c) and by 
symmetry a;+(c) play the role of bifurcation parameters for the stochastic 



5 



FitzHugh-Nagumo system fll.3l) -(1.4). Indeed for a in the neighborhood of 
x_(c) but smaller the limit of (X^, Yt) is a unique equilibrium point, whereas 
for a in the neighborhood of X-{c) but greater it is the graph of a periodic 
function. These bifurcation parameters are different from those of the de- 
terministic system (11. ip . This theorem is a new result w.r.t. [6J. It is made 
possible by the freedom on parameter a. The new limit cycle '^'^{t)) 
is defined in the same way as in [B] Theorem 1 Part 3., provided we take into 
account the presence of a in our system. 

Other regimes are considered in the work of Berglund and Gentz (cf. |JJ and 
references therein). 

Theorem 2.2 Let c > S . Consider x*_{y*) and x\{y*) defined in proposi- 
tion 3.2 and definition 3.3; then for all y €]/(ao), /(ai)[ and for all h > 
there exists t{y, h) such that for all A > t{y,h), 

1. Ifae]x*_{y*),xl{y*)[ 

lim*P(,,^)( sup \Yt-y*\>h) = (2.5) 
[t(.y)A] 

2. If a < x*_{y*) or a> x\{y*), 

lim*P(,,,)( sup {\Xt -a\ + \Yt - /(a)|) > h) = Q (2.6) 

Remark: 

Case 1 may be considered as a degenerate version of the limit cycle of case 1 
theorem 2.1. In fact y* is a fixed point but Xt oscillates between x*_{y*) and 
x\{y*). 

3 Exit Time, Main State and Metastable State 

3.1 Basic results on Large Deviations Theory 

Because of the slow-fast property of FitzHugh-Nagumo systems, the slow 
variable Yt of system (1.3)-(1.4) may be in a first approximation frozen at 
the value y which leads us to the study of the family of one dimensional 
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dynamical systems indexed by parameter which plays a basic role in the 
study of FitzHugh-Nagumo systems ([III])-(1.2) and (OD-fl.4): 

dxt = {-y + f{xi))dt, xl = x (3.1) 
So we are led to consider the real valued deterministic system 

dxt = b{xt)dt, Xq = X (3.2) 
and its perturbation by a brownian motion 

dxt = b{xt)dt + VidWt, xq = x (3.3) 

Let us briefly recall some results from [5J. The process [xt] describes the 
movement of a particle on the real line submitted to the force field h{x) and 
to a stationnary Gaussian noise of amplitude v^l. When e ^ 0, (x^) converges 
to the solution {xt) of fl3.2p : 

Vr7>0 VT>0 limP(sup|xt-Xt| > r/) =0 (3.4) 

[0,T] 

However because of diffusion due to the presence of noise, some trajectories of 
the process {xt) may present large deviations from those of the deterministic 
system (xj). Such deviations are measured by means of the action functional 
S'^l{ip) independant of e and defined by 

StU'^) = o / \'fu- h{ipu)\^du (3.5) 

when ip is absolutely continuous, by S'^p^^ip) = +oo otherwise 

Theorem 3.1 ( l^, Lemma 2.1, Chap. 4) 
Let 1] > 0. Then 

Pa;(sup \xt - xt\ >r])< exp(-^[inf S'^(93) + o(l)]) (3.6) 

[o,T] e ^ 

when e ^ and where A := {(p; (p^ = x, sup^j.] \ipt — Xt\ > rj}. 

Large deviations theory also provides estimates on the first exit time of (xj) 
from a domain (cf. [5J, Theorem 4.2, Chap. 4). Domains of interest are 
basins of attraction of stable equilibrium points of (13. 2p . The key tools are 
quasipotentials. 
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Definition 3.1 The quasi potential of the deterministic system liS. 2\] w.r.t. 

a point X (also called transition rate) is defined as the function 

u ^ V{u) := infiS^li^y, < Ti < T2, <^(Ti) = x, (^(Ts) = u} (3.7) 

Proposition 3.1 The quasi potential of ( 15*.^)) w.r.t. x coincides with the 
function 

ru 

u ^ V{u) = -2 /_ b{r)dr (3.8) 

Remark: 

The above statement holds since (13. 2p is one dimensional. It also holds in 
the multidimensional case when the drift h of (13.21) is a gradient. 

Theorem 3.2 Let x* he a stable equilibrium point of liS. 2\] such that b{r) < 
for all r > X* , and b{r) > for all r < x* . Let D be the basin of attraction 
of X* and f denote the first exit time of x from D. Let us assume that 
D =]«!, with V{ai) < V{a2). Then for all x & D and h > 0, 

limP^(i^ = ai) = 1 (3.9) 

V(Qi)-h V{a-,) + h 

limP^(e^^ < f < e^^) = 1 (3.10) 

Remark: 

With the notations of Theorem 13.21 with great probability when e — 0, the 
behaviour of the process (xt) on an interval [0,T(e)], in particular whether 
the process has jumped out of the basin of attraction D or not, depends on 
the value of elogT(e) compared to V{ai). 



3.2 Two Families of Quasipotentials. 

Let us now apply these results to the following family of one dimensional 
dynamical systems indexed by parameter y introduced in the preceeding 
subsection. 

dx! = i-y + f{x',))dt, xl = x (3.11) 

Proposition 3.2 Let y G]/(cto); /(c^i) [ with and a\ the two points where 
f vanishes. 

(i) The set {x G R; f{x) = y} consists of three points x*_{y) < XQ^y) < x^(y) 
( the equilibrium points of l{3.11\) ) each being a continuous function of y with 
bounded first and second derivatives, 
(a) The two points x*j^{y) are stable, x^i^y) is unstable. 
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Proof of Proposition 13.21 

Left to the reader; we refer to figure 1-section 2 



Definition 3.2 Let us define the two functions V± on ]/(ao), /(cti) [ as fol- 
lows: 



From Proposition 13. II we see that V±{y) is the value at xliij) of the quasipo- 
tential of (13.21) w.r.t. x^(?/). Both functions V±{y) are strictly monotone: 
V- is strictly increasing, Vj^ is strictly decreasing. Therefore their graphs 
restricted to ]/(ao), /(ai)[ intersect at a unique point. 

Definition 3.3 We denote by {y*, S) the intersection point of the graphs of 
K_ and V+. Let Si := {y > y*} and £2 '■= {y < y*}- 

Definition 3.4 For c g]0, S[ we denote by y±{c) the points o/]/(ao), /(a.i)[ 
which satisfy y^{c) < y* < y+{c) and V_(y_(c)) = c = V+(y+(c)). Let us 
also define X-{c) := x*_{y^{c)) and x+(c) := x\{y+{cj) (cf. figure 1). 



By definition V-{y*) = V+{y*) = S. For f{x) = Ax - x^, y* = and 3 = 4. 
For any function U satisfying y — f{u) = —dxU/2 the following identities 



V^y) = U{x*{y)) - UixUy)), V+iy) = U{x*{y)) - U{x\{y)) (3.13) 



In our case such a function f/ is a polynomial of degree 4 which admits 
Xq(?/) as relative maximum and x*^{y) as relative minima. The graph of the 
function U has two wells with respective bottoms at x^(y) and one top at 
xliy). Identities (13.131) express that V±(?/) are the respective depths of these 
wells. Therefore, on £1 the well with bottom x*_{y) is the deepest one while 
it is the contrary on £2. 

Second, the portion of the curve z = f(x) connecting {x*_(y), y) to (xQ(y), y) 
is situated below the horizontal line Ly := {(x, z); z = y}; thus the positive 
quantity ^V_{y) measures the surface of the area limited by the curve z = 
f{x) and the segment of Ly on which x G [x*_{y),XQ{y)]. In the same way 
^^+{y) measures the surface of the area limited by the curve z = f{x) and 




(3.12) 



Remark: 



hold: 
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the segment of Ly on which x G [x^iy) , x*^{ii)\ but in this case the portion of 
the curve is above the hne segment. 

Moreover as we will see in the following section V± is connected to exit times 
of diffusions 

d~Zl = {-y + fCZ'!))dt + VldWt, Zl = X (3.14) 
from the basins of attraction of x*j^{y) (cf. Theorem 13.21) . 

3.3 Exit Times, Main state and Metastable States 

We refer the reader to [7] for the present section. Let us recall the fundamen- 
tal difference between the two parameters e and 5. Parameter 5 is already 
present in the deterministic system fll.ip -(1.2) where it measures the differ- 
ence between the time scale of the slow variable yt and the time scale of the 
fast variable Xt- In particular after the time change s := 5t, the trajectory 
{xs.ys) ■= {xst.yst) satisfies 

Xs = -ys + f{xs) (3.15) 
y^ = S{xs-a) (3.16) 

Since 6 is small the component jjs may be considered as constant equal to yo. 
Let us define to be the family of solutions of 

SdZf = {-y + f{Zl))dt + ^tdWu Zl = X. (3.17) 

We will also consider the time change {Zf) of (Zf) under s := 5t 

dZl = {-y + f{Zy))ds + VidWs, A = x (3.18) 

with e = |. (Zf) is the stochastic perturbation of (13.111) . 

In the two following propositions we give crucial estimates on exit time of 

the respective solutions of (3.17) and (3.18). 

Proposition 3.3 Let ff (resp. ff ) denote the exit time of Z^ from Df (resp. 
D2) which is the basin of attraction ofx*_{y) (resp. x*j_{y)). Let us recall that 
D\ =] — oo,XQ(y)[ (resp. =]3^o(2/)' +^[/ From Theorem \3.^ identity 
i\3.10\) . for X G Df and h > we obtain 
Wx g] - oo,x*{y)[,yh > 

limP.(exp(^^M^) < ff < e,p( ^-(^)) + ^ )) = 1 (3.19) 
e^o e e 
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An analogous result holds for x E D2 by replacing rf by rf and V^{y) by 

Vx e]x*Q{y),+oo[,Wh > 0; 

limP.(exp(^^tM^) < fv < exp(^^M)±Z^)) = 1 (3.2O) 

e^o e e 

Proposition 3.4 Let rf (resp. T2) denote the first exit time of from D\ 
(resp. 02). The law of rf is the same as the law of Sff fori = 1,2. Let us 
assume that e and 6 go to in such a way that 

Illegal ^cG]0,+oo[. (3.21) 


In this case we obtain: 
\/x g] - 00, xl{y)[,\/h > 

limP,(5^"'(^-^-(^)+'^) < rf < = 1 (3.22) 

An analogous result holds for x G D2 by replacing rf by rf and V^{y) by 

V+{y): 

Vx g]x;5(?/),+oo[,v/i > 

limP,(5'="'('=-^+(^)+'') < r| < S'^-'i'^-y+iy)-h)) = 1 (3.23) 

These two properties enable us to introduce some remarks about the 
main theorems stated in the previous section. These remarks are hnked to 
the notions of mainstate and metastable state introduced in [5]. 
In a general framework ( cf. [7]) the main state is the point towards which 
the cost of moving, or the transition rate, is minimum. It is not always 
unique: for instance in our bistable case there are two main states when 
y = y*. Before reaching the main state, the process may reach metastable 
ones accessible for shorter time lenghts. 

Main states may be considered as stable states, whereas metastable states 
are only stable in some time scale. 

To study transitions of fl3.18p between the two basins of attraction dur- 
ing [0,T(e)], the relevant quantity to consider is elogT(e) which we must 
compare to the transition rates V±{y) given by fl3.12l) . For the time scale 
T(e) = e^; elogT(e) = c. So we must compare c to the transition rates 
V±{y). Actually this amounts to compare first c to 5* defined in definition 
3.3. We refer again to figure 1. 

More precisely we can state using definition 3.3 and 3.4: 
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Proposition 3.5 ; 1. When c > S the main state of is equal to x^{y) 

(resp. x*_{y)) if y < y* (resp. y > y*). 

When y = y* the two points a;^(y*) are both main state. 

2. When c < S for y E]y-{c),y+{c)[ and x & Df (resp. x G D2) the 

metastahle state of is equal to x*_{y) (resp. x*^[y)). 

Proof of Proposition 13.51 

Direct consequence of the estimation of the time of exit given in proposition 
3.4. 

Remark: 

1. When c> S the time interval [0,T(e)] is long enough so that the process 
Zy reaches with great probabihty a small neighborhood of its main states. 
And as we can find an open interval / containing y* such that: 

\JyeI c>max(Vl(y),l^+(|/)). (3.24) 

following proposition both exit time rf and rf tends to in probability so 
we can only expect results on the slow component and the result depends on 
whether the boundary between the two main states is attracting or notas is 
shown in theorem(2.2). 

2. On the countrary when c < S the time interval [0, T(e)] is too short and the 
process only reaches with great probability a neighborhood of a metastable 
state. In this case we have y-{c) < y+{c) and if we consider the interval 
]?/„(c), t/+(c) [ when c < V^iy) < V^iy) and x G D\ (resp. c < V^iy) < 
V ^ — {y) and x G -Df)) ^lU y G]?/_(c),y+(c)[ exit time rf (resp. )tends 
to infinity. So the process Z^ remain in a neighborhood of one metastable 
state and switch to the other one as soon as Yt gets out of ]y_(c), |/+(c)[; this 
give rise to a limit cycle as is shown in theorem (2.1). 



4 Proof of the Main Theorems 

Definition 4.1 For a e]x-{c).,Xj^{c)[ and y G]y-(c), y+(c) [ as we can check 
on figure 1 we have: 

x*^{z) < a;„(c) <a < x+{c) < x\{z) (4.1) 

Then for c e]0,S[, a g]x_(c), x+(c)[ and y g]?/_(c), ?/+(c)[, the following 
definitions of time duration makes sense 

T^(c) = / and T^ic) = / r (4.2) 

' Jy-(c) xUy)-a ' JyAc) \xt{y)-a\ ^ ' 
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where y±{c) and x±{c) have been introduced in Proposition 3.2 and Definition 

Z3 



On the interval [0, H-Tg^jwe can introduce the solution of the following 
ordinary differential equations used to define the limit cycle in Theorem 

Definition 4.2 For c &]0,S[, a g]x_(c), x+(c)[ and y g]?/_(c), y+(c) [ we 
define as the continuous periodic function with period T"-{c) = (c) + 
T2"(c) satisfying \E'c(0) = y-{c) and the ode 

Kit) = ximt))-a, tG[0,Tf(c)[ (4.3) 
vi>«(t) = x*_mt))-a, te[T,^{c),T,^{c) + T,^{c)[ (4.4) 

Fort ^ {kT'\c),kT''{c) +T^{c),k E Z}, we denote by ^'^(t) the derivative 
of'^"^ att. 

Existence and regularity of solutions to odes (4.3)-(4.4) follow from the 
regularity of the functions x±.. Under the assumptions we recall that we 
have (see figure 1) 

x*_{z) < x-{c) < a < x+(c) < x*^{z) (4.5) 

so \E'c(t) increases from y-{c) to y+{c) in a duration of time equal to T{'(c) 
and decreases from y+{c) to y_(c) in a duration of time equal to T2 (c). The 
periodic function \E'" is obtained by continuously sticking together at y±{c) 
the solutions y°i_ of the following odes 

Vlit) = ximt)) - a, r±{0)=y (4.6) 

Remark: 

The function \E'" is continuous on R differentiable with continuous derivatives 
except at points in the set {kT^ic), kT^ic) + Tf (c), k e Z}. 
For t e [0,rf(c)[ the point (x:;(^;!(t)), ^;?(t)) belongs to the right stable 
branch; for t E [Tf{c), Tf{c) + T^lc)[ the point (xl(^;?(t)), ^:?(t)) belongs to 
the left stable branch. 

Proof of Theorem [H 

1. Case a G]a;_(c), x+(c)[: 

a) Let y < y-{c) and x G Di. Then lim*P^(r|' = 0) = 1 since V^{y) < 
y_(?/_(c)) = c (remember that is strictly increasing). Therefore the 
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process (Xt) leaves Df instantaneously and is attracted to a neighborhood 
of x'^ly) > a;+(c). The point {Xt,Yt) remains in a neighbourhood of the 
branch v = f{u) containing {x^{y),y) and Yt increases as long as Xt > x~^{c) 
since dVt = {Xt — a)dt. However identity (3.23) implies that lim*P3;(r| = 
+00) = 1 for z < y+(c) (resp. lim*P2,(r| = 0) = 1 for 2; > y+{c) ). 
Therefore the point {Xt, Yt) is instantaneously attracted to a neighborhood 
of (xl(y+(c)), y+{c)) after Yt has crossed y+{c) since the speed of Yt is strictly 
positive in a neighbourhood of (x+(c), ?/+(c)). The argument is then the same 
as before. We detail it for completeness. Since x*_{y+{c)) < x_(c) < a, the 
second coordinate Yt decreases as long as Xt — a < 0. However identity (3.22) 
implies that lim*Pj.(rf = +00) = 1 for z > y-{c) (resp. lim*Pa,(rf = 0) = 1 
for z < y_{c)). Therefore the point {Xt,Yt) is instantaneously captured by 
(x^(?/_(c)), y_(c)) after Yt has crossed y-{c). Hence {X,Y) converges in 
probability to a limit cycle of period T"(c) = (c) + (c). 
b) Let y > y-{c) and x E D\. By identity 3.22 lim*Pj;(rf = +00) = 1, the 
fast process {Xt) is attracted to x*_{y). After this first phase, one can apply 
the same argument as in a). 

2. Case a < X-{c). As already mentioned in section 3, when c < S, the 
relevant states are the metastable ones. The assumption a < x_(c) implies 
f{a) > y^{c). Let V be a small neighborhood V of {a,f{a)). Let us first 
assume that V is so small that all {u,v) G V satisfies v > y-{c) and ac- 
cordingly u < x_(c). Let {Xt,Yt) start from {u,v). In particular u E D\ 
and tI = +00 since V-{v) > c. Then {Yt) will evolve like the solution of 
Vt = x*_{vt) — a, vo = V for which f{a) is a stable equilibrium point. Let us 
now assume that there are points {u,v) G V such that v < y-{c). Then {Xt) 
instantaneously jumps to the right stable branch; {Yt) becomes close to the 
solution of Vt = x^{vt) — a and therefore increases until it reaches y+{c). At 
that time it jumps to the left stable branch and we are back to the previous 
argument since then {Yt) becomes close to the solution of Vt = x*_{vt) — a. 
Let us notice that for the slow variable Yt we get a result using the uniform 
norm by large deviation estimates but for the quick variable Xt the result is 
formulated in norm thanks to equation (1.4). 

Proof of Theorem 12.21 When c > S", as pointed out in subsection (3) 
has time to reach its main state so the evolution of Yt is close to the solution 
of 

Vt = x*^{vt) — a provided Vt < y* 



14 



Vt = x*_{vt) — a provided Vt > y* 

and it depends on whether the boundary {y*} between Si := {y > y*} 
and £2 := {y < y*} is attractive for this system or not(cf. [7]). If a G 
]x*_{y*),x*^_{y*)[ the boundary is attractive. This is not the case when a < 
x*^_{y* nor when a > x*^_{y*). 

1. Let x*_{y) < a < x\{y*). The point y* is an attracting boundary since 
x*_{y*) - a < < x\{y*) - a. 

Assume for example that Yq = y < y* that is ?/ G 82- The evolution of (Yt) is 
close to the solution of Yt = x'!\_{Yt) — a. Therefore since x'\_{y) > a, process 
(Yf) starts increasing until it reaches y*. After it has reached y*, the evolution 

of (Yt) is close to the solution of Yq = y* , Yt = B(Yt) where -B is a vector 
field tangent to boundary (cf. [7j Theorem 1) which is the 0-dimensional 
manifold {y*}', therefore B is zero and Yt = y*. The time t{y, h) is the time 
necessary to reach a small ball B{y*, h). 

\i y & £1 the evolution of {Yt) is close to the solution of 3^t = x*_{yt) — a. In 
this case process (Yt) starts decreasing since < a until it reaches y*. 

The argument and the conclusion are then identical to the preceeding case. 

2. Let a > a;^(y). The point y* is not an attracting boundary since x*_{y*) — 
a =< and x*^{y*) - a < 0. 

Let us assume for example that Yq = y E 82- The evolution of (Yt) is close 

to the solution of Yt = x!j_(yt) — cl- The value /(a) is a stable equilibrium 
for (Yt). If x'^iy) > CL, process (Yt) starts increasing until it reaches /(a) to 
which it is attracted. If now a;^j_(y) < a the argument and the conclusion are 
the same except that (Yt) starts decreasing, li y E Si the evolution of (Yt) is 

close to the solution of yt = x*_{yt) — 0,. Since x*_{y*) — a < process(yt) 
starts decreasing until it crosses y* (which is not attractive) towards 82- Then 
we are back to the previous case. 

3. The case a < x*_{y) is treated analogously. 
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